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ABSTRACT 

We study how rotation affects observable amplitudes of high-order g- and mixed r/g-modes and examine 
prospects for their detection and identification. Our formalism, which is described in some detail, relies on a 
nonadiabatic generalization of the traditional approximation. Numerical results are presented for a number 
of unstable modes in a model of SPB star, at rotation rates up to 250 km/s. It is shown that rotation has a 
large effect on mode visibility in light and in mean radial velocity variations. In most cases, fast rotation 
impairs mode detectability of g-modes in light variation, as Townsend (2003b) has already noted, but it helps 
detection in radial velocity variation. The mixed modes, which exist only at sufficiently fast rotation, are also 
more easily seen in radial velocity. The amplitude ratios and phase differences are strongly dependent on the 
aspect, the rotational velocity and on the mode. The latter dependence is essential for mode identification. 

Stars: oscillations - Stars: emission-line, Be — Stars: rotation 

1 Introduction 

Variability with frequencies comparable to rotation frequency has been found in a num- 
ber of hot (most often Be) stars. Whether such variability is caused by slow modes has 
been debated for some time (see e.g., Baade 1982, Balona 1985). However, recent data 
from the MOST on £ Oph (Walker et al. 2005a), HD 163868 (Walker et al. 2005b), 
and p CMi (Saio et al. 2007) revealed rich frequency spectra which may be understood 
only in terms of oscillation mode excitation. There are also frequency spectra obtained 
from ground-based observations, such as of fj Eri (Jerzykiewicz et al. 2005), which 
certainly cannot be explained in terms of a rotational modulation. A potential for using 
the abundant frequency data as constraints on stellar models exists but the prerequisite 
is mode identification and it has not been done so far. 

Most frequently used method of mode identification employs the data on mode 
amplitudes and phases of light variability in various bands and in radial velocity. Am- 
plitude ratios depend on angular dependence of surface distortion, which in rotating 
stars is no longer described by a single spherical harmonic. The departure becomes 
large once the angular velocity of rotation is comparable to oscillation frequency. This 
is not an unusual situation for cooler B stars in the main sequence band. Take a model 
of 6 M star in the mid of its main sequence evolution. Ignoring effects of rotation, we 
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find that all dipole modes with period between 1 .8 d to 2.8 d are unstable. Rotation pe- 
riods in this range correspond to equatorial velocities between 100 km/s and 150 km/s, 
which is not high for such a star. 

For low frequency modes the surface dependence may be approximately described 
in terms of the Hough functions (see e.g., Lee and Saio 1997, Bildsten et al. 1996, 
Townsend 2003a). This approximation, called traditional, was used by Townsend 
(2003b), who calculated observable light amplitudes and used the results to address the 
problem of mode identification for low frequencies in rotating stars using multicolor 
data. It is our experience (e.g., Daszyhska-Daszkiewicz, Dziembowski and Pamyat- 
nykh 2005), however, that in the case of B-type pulsators, it is very important to com- 
bine photometric and radial velocity data for a unique discrimination of excited modes 
and constraining stellar parameters. Thus, this work, which may be regarded as an 
extension of Townsend's paper, focuses on calculation of radial velocity amplitudes. 
We adopt an uniform approach in our calculation of all disk-averaged amplitudes and 
it is different from that used by him. Furthermore, in addition to g-modes, we include 
r-modes, which become propagative in stellar radiative envelopes once rotation is fast 
enough (see Savonije 2005, Townsend 2005b, who uses the term mixed gravity-Rossby 
modes, Lee 2006). 

In Section 2, after specifying assumptions adopted in our calculations, we summa- 
rize formulae for angular dependence of velocity and atmospheric parameters. Expres- 
sions for the light and disk-averaged radial velocity variations are given in Sections 3 
and 4, respectively. Sections 5, 6 and 7 present numerical results for selected modes 
in one representative stellar model considering range of rotation rates but ignoring ef- 
fects of changes of centrifugal force on model structure. Unstable mode properties 
are briefly described in Section 5. In Section 6, upon adopting an arbitrary normal- 
ization of linear eigenfunctions, we calculate observable amplitudes of various modes, 
which may be excited and detected. Prospects of mode identification are discussed 
in Section 7. Examples of diagnostic diagrams employing amplitude ratios and phase 
differences are shown there. 



2 Photospheric Parameter Variations and Velocity Field 
in the Traditional Approximation 



In our study we adopt the standard approximations, which are the linear nonadiabatic 
theory of stellar oscillation and static plane-parallel atmosphere models. These ap- 
proximations are well justified in our application. Modes detected in slowly pulsating 
B-type stars have indeed very low amplitudes, the oscillation periods are much longer 
than the thermal time scale in the atmosphere, and vertical variations of mode am- 
plitude are small across the whole atmosphere. The constant kinematic acceleration 
is easily included. Like Townsend (2003b) we adopt the traditional approximation, 
which allows to separate latitudinal and radial dependencies of the pulsational ampli- 
tudes. We essentially follow his formalism, expect that we choose the azimuthal and 
temporal dependence as Z = exp[i(mcp — COf )], which implies m > for prograde modes 
and m < for retrograde modes. 

The displacement vector at the surface may be written as follows 




(1) 



where e is an arbitrary, but small, complex constant, and 



GM 



CO 2 /? 3 ' 
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The three functions ©, 0, describe the latitudinal dependence of solutions and are 
the Hough functions, which are obtained as solutions of the Laplace's tidal equations 

(©+/jmju)0=(j 2 ju 2 -1)0, (2a) 
(®-m S( u)0=[A,(l- A / 2 )-m 2 ]0, (2b) 

where s = 2Q./& is called the spin parameter, /j = cos8 and <D = (1 — A* 2 )^"- The 

equations together with boundary conditions at /j = and /a = 1 define the eigenvalue 
problem on X. The third function is given by 

® = -m® + sn®. (3) 

Bildsten et al. (1996), Lee and Saio (1997), and Townsend (2003a) discussed in great 
detail the X{s) dependence and asymptotic properties of the Hough functions. Here, 
we recall only the essentials. 

For specified m and s^O, there are branches with X^>£(£+ 1 ) and 0(8) -^pj m ' (ex- 
cept for normalization). These branches correspond to g-modes distorted by rotation. 
For prograde sectorial modes (m=£), X slowly decreases with s. For all other g-modes, 
the function X(s) is increasing quite rapidly. We will identify g-mode branches by the 
£ value at s = 0. Thus, like in the case of no rotation, we will use (£, m) values as the 
angular quantum numbers and refer to £ as the mode degree. However, now specifi- 
cation of the angular dependence requires also the value of s. The symmetry about 
the equatorial plane is determined by the parity of £+ \m\. If it is even, the functions 
0(8) and 0(8) are symmetrical and 0(6) is antisymmetrical. The opposite is true if 
£+\m\ is odd. The branches for which X — > — °° at s — > correspond to r-modes. For 
each m < 0, there is one branch crossing zero at s= \m\ + 1. If X > 0, the associated 
modes become propagatory in the radiative regions, they may be excited and visible in 
the light variations. However, following Lee (2006), we will still call them r-modes. 
The functions 0(8) and 0(8) are antisymmetrical and 0(8) is symmetrical about the 
equator for these r-modes. 

Upon replacing £(£+!) with X(s), the nonadiabatic mode properties may be calcu- 
lated with a reasonable accuracy using the same code as for non-rotating stars. This is 
so because for the mode of our interest, horizontal flux losses, which are not correctly 
described, are small. The latitudinal dependence of perturbed thermodynamical pa- 
rameters is then given by 0(8). Thus, the bolometric flux perturbation may be written 
as 

^= £/ ©Z, (4) 

^bol 

where / is a complex quantity determined by solution of linear nonadiabatic equations 
and it depends on the stellar model and the mode parameters m, X and co. 

For evaluation of perturbed monochromatic fluxes, we also need the perturbed grav- 
ity, which, as follows from Eq. (1), is given by 

^ = -e(2 + O3- 1 )0Z. (5) 
8 

The effect of gravity perturbation plays a relatively small role in light variability caused 
by slow modes but it is easy to include. More important is perturbation of star shape 
leading to changes in the projected surface element and the limb-darkening. For both 
we need the normal to stellar surface which, as follows from Eq. (1), is given by 

5 „, = - E V„ (9Z ) = - e (0»,^)z. (6) 
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The change of the directed element of the surface is 

SdS _ / 3© _v 

~ds~~ e ( ' W sine 

where dS = /? 2 d/j dep. 

From Eq. (1), we also obtain the perturbed pulsation velocity field as seen from an 
inertial system (90 = <P + ^0 

Use of the Lagrangian pulsational velocity is adequate for representing velocity at the 
photospheric layer because % is nearly constant across the outer layers for the slow 
modes considered by us. 



3 Light Variation 



3.1 Semi- Analytical Formula 

The most straightforward extension of the expression for the light variation to the case 
of rotating stars is through the expansion of the Hough function into the truncated 
series of the associated Legendre functions. This was the way Townsend (2003b) 
derived his expression. We write the equivalent expression in the form, which is a 
straightforward generalization of our formula (Daszyhska-Daszkiewicz, Dziembowski 
and Pamyatnykh 2003) derived for the case of modes described by single spherical 
harmonic. Now the complex amplitude of the light variation in the x passband may be 
expressed as 

(9) 



where 



and 



7=1 

H +2(7-1) 
|m|+2(7-l) + l 

©? = -1.086fc 



<D 



-1.086&J 



(2 + £)(!-£)- (2 + (S 



even — parity modes 
odd — parity modes 

l 9log(j,|^|) 

4 3iogr e ff 

_ n aiog(^|) 



dlogg 



b*e = f fhWpPefrW. 



(10) 



where h x is the limb-darkening law, adopted in the nonlinear form (Claret 2000). With 
this form Townsend obtained his analytical expressions for M. The quantities 7? (s), 
which have to be calculated numerically, are the expansion coefficients of the © func- 
tion into the series of the Legendre functions. That is 



0(6) = 1^.(^(6). 



(11) 



The expression given by Eq. (9) is quite revealing. At specified e, the light ampli- 
tudes of low degree modes must decrease with s because of increasing role of higher 
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order terms which suffer more from disk-averaging. The higher order terms lead to the 
aspect-dependence of the amplitude ratios. 

Unfortunately, we could not find a corresponding semi-analytical expression for 
the radial velocity and this is why we decided to rely on two-dimensional numerical 
integration over the visible hemisphere. We used Eq. (9) to check the accuracy. 

3.2 Numerical Approach 

The total flux in the x passband toward the observer is given by 

L x = J JAnobs-dS (12) 
s 

where integration is carried over visible part of star surface, S, and n b s is the unit 
vector toward observer. 

In the spherical coordinate system with the polar axis parallel to the rotation axis, 
we have 

n D bs= K,0e,O(p) (13) 

where 

o r =Ji — cos / cos + sin / sin 8 cos (tp — <po), 

09 = — cos /sin 8 + sin /cos 8 cos (cp — (po), 

Oq, = — sin/sin(cp — (po)- 

The observer's angular coordinates are (Z,cpo)- The first order perturbation of the total 
flux is given by 

8L X = J [(8 J A + TMx) dSe r + T x h x 8dS] • n obs . (14) 
s 

Assuming an equilibrium atmosphere, we have from Eqs. (4) and (5) 

0$ 



5j 
where 



:/-o^(2+ra- 1 ) 



©z 



«*r = |^ and t£ = ^ 



d\ogT eS 8 3logg 

are the derivatives determined numerically from grids of stellar atmosphere models. 
For perturbed limb darkening, we take into account perturbations of the coefficients 
on r e ff and g though they are only secondary contributors to light variation. More 
important contribution arises from 8n y (see Eq. 6). With all terms included, we obtain 

(\dh x f dh x , ,J dh x fd& im© 

The derivatives of Claret's h x are given in Appendix Al. The expression for 
perturbed surface element is given in Eq. (7). The integration is carried over the unper- 
turbed visible hemisphere. Within the linear approximation, the integration boundary 
is unchanged. Note also that the horizontal component of the displacement does not 
enter the expression. The domains of integration over 8 and (p are shown in Fig. 1 and 
Fig. 2, respectively. From Figs. 1 and 2 it follows that the ranges are: 



0<8<! + /, 
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( obs 




- (71 — a) < cp - cpo < {n - a) 

where a = arccos[cot8cot/]. For each 8, we carry the integration over the azimuthal 
angle, x P = (p— cpo, from — p to +p, where p = 7C— a = 7C— arccos[cot8cot/] (see Fig. 1 
and Fig. 2). We use identities 



r 

/ 



2jG( l P)cosm l Pd l P if G(»P) is even 
o 

P 

2ij'G( v I')sinm l I'd<f l I' if G(¥) is odd 



The final expression for the total light variation can be written in the following form 



8L X 



a; 



®i + ®3 )/ + 2«i + « 2 -(2 + 03 1 )(aJ«i + « 4 



U5) 



and Zo = exp[i(m(po — CM)]. In Appendix A2 we give explicit expressions for the two- 
dimensional integrals, <B , which depend on two angular numbers, spin and the aspect. 



7 



The integrals take into account changes in the limb-darkening resulting from the change 
of the normal (Eq. 6) as well as the change due to perturbation of the local temperature 
(Eq. 4) and gravity (Eq. 5). The two latter are given through derivatives of h x with 
respect to logoff and logg. There are many terms contributing to 8L X . However, in our 
applications two are far dominant: the one resulting from temperature perturbation, 
which is proportional to OCySi, and the other resulting from the surface distortion, 
which is proportional to $2- 



4 Radial Velocity Variation 

Adopting the standard sign convention, we write the radial velocity averaged over the 
stellar disk in the following form: 



(V rad > = 



f(v-n obs )j /z x n obs -dS 

Jj> 

/ f h x n obs • dS 

So 



where for the total velocity field we use 

v = 8v + i2/?sin0e<p. 
The pulsational component, 8v, as results from Eqs. (1) and (8), is given by 



(16) 



(17) 



8v = 



0--C5© 

2 



\ 



-iC3 



03 











S COS0Q 



sin 9 2 sin 8 

S COS0 



sin0 2 sin0 







+ i-0sin9 



(18) 



The contribution of rotation to the mean radial velocity arises from the same pulsational 
changes of photospheric parameters which cause luminosity change and may be calcu- 
lated in the same way as outlined in Section 3. Clearly, there is a nonzero contribution 
only for non-axisymmetric modes. 

We write our final expression for the perturbed radial velocity in the following 
form, 

8(y rad ) = icoe/?(c P uis + Got)Zo, (19) 



with 
and 



Cn 



ssmi 



Cpuls = (®5 - l^l) + 05 (®6 - |«s) 

^af!89 + Sii )/+2®9 + «io- (2 + C5- 1 ) (a^« 9 + «12) 



Explicit expressions for the <£ coefficients in terms of the three Hough functions are 
given in Appendix A3. 

If s w 1 and m^O, the two contributions are of the same order. Similarly to the 
case of the luminosity, the dominant contributions to C m t are the ones resulting from 
temperature perturbation, which here is proportional to aj'Sg, and the other resulting 
from the surface distortion, which is proportional to ®iq. 
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5 Unstable Modes in a 6 M Main Sequence Star 

To illustrate how visibility of various modes depends on equatorial velocity of rota- 
tion we choose a model of a 6 M Population I star in the mid of its main sequence 
evolution. Parameters of the model are given in Table 1 . The model is spherically sym- 
metric, which is consistent with our use of the traditional approximation, but includes 
averaged effects of centrifugal force corresponding to uniform rotation with equatorial 
velocity of 250 km/s. The mean effects of centrifugal force are still reasonably small 
and therefore we used the same model to study effects of the Coriolis force at lower 
equatorial velocities. 

There are many low frequency modes that are unstable in our selected star model. 
These are mainly g-modes. However, at sufficiently high rotation rate there are also 
certain r-modes, which become propagatory in radiative envelope and may become un- 
stable. In all cases, the instability is caused by the K-mechanism acting in the metal 
opacity bump layer. For g-modes, the angular order, £, is defined at the limit of zero 
rotation where the angular dependence of amplitude is described by individual spher- 
ical harmonics. The actual mode geometry is determined by the (£,m) numbers and 
the spin parameter, s. Since at each azimuthal order m < there is only one r-mode for 
which X changes sign, we will identify r-modes by the m value alone. Naturally, we 
focus attention on modes suffering least reduction of observable amplitude caused by 
disk-averaging. Therefore, we consider g-modes with £ < 2 and r-modes with m=—\ 
and —2. 




e e 

Fig. 3. The Hough functions for the selected g-modes with (. = 1 and r-modes with m = — 1 (see Table 1 for 
more data on the modes). 

Typically, at each degree and azimuthal order, we find instability extending over 
many (up to 40) radial orders. For analysis of visibility, we selected the mode charac- 
terized by the highest normalized growth rate, T), which varies between —1 and 1. The 
important parameters of the selected modes at adopted equatorial velocities are listed 
in Table 1. For the r-modes we write z"instead of £ as the first entry. The effect of 
rotation on mode surface geometry at specified V rot depends on s which determines X. 
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Table 1 

Most unstable low degree modes in the B star model with M = 6.0 M©, logoff = 4.205, logL/L = 3.204, 
for the four values of rotational velocity 



I 


m 


Vrot 


spin 


X 


CD 


/ 


Vobs [c/d] 


Vstar [c/d] 


Tl 


1 








0.00 


2.00 


17.76 


(9.99, 11.55) 


0.4220 


0.4220 


0.110 


1 





50 


0.86 


2.33 


16.35 


(10.56, 12.18) 


0.4396 


0.4396 


0.129 


1 





150 


1.97 


4.62 


9.48 


( 8.55, 13.02) 


0.5772 


0.5772 


0.202 


1 





250 


2.65 


7.61 


6.16 


( 5.82, 12.46) 


0.7162 


0.7162 


0.238 


1 


+1 





0.00 


2.00 


17.76 


(9.99, 11.55) 


0.4220 


0.4220 


0.110 


1 


+1 


50 


1.02 


1.45 


22.70 


(10.38, 10.87) 


0.5627 


0.3731 


0.063 


1 


+1 


150 


3.27 


1.17 


26.11 


(10.00, 10.19) 


0.9167 


0.3478 


0.031 


1 


+ 1 


250 


5.63 


1.09 


27.82 


(10.33, 10.14) 


1.2851 


0.3369 


0.021 


1 


_1 





0.00 


2.00 


17.76 


(9.99, 11.55) 


0.4220 


0.4220 


0.110 


1 


_1 


50 


0.72 


3.60 


11.37 


(8.98, 12.61) 


0.3375 


0.5271 


0.179 


1 


-1 


150 


1.32 


11.80 


4.25 


(3.42, 11.27) 


0.2930 


0.8619 


0.258 


1 




250 


1.71 


22.21 


2.57 


(0.42, 8.50) 


0.0870 


1.1097 


0.265 


2 








0.00 


6.00 


7.35 


(6.46, 12.38) 


0.6569 


0.6569 


0.225 


2 





50 


0.54 


6.84 


6.42 


(5.41, 11.92) 


0.7015 


0.7015 


0.234 


2 





150 


1.24 


13.37 


3.76 


(2.43, 10.46) 


0.9170 


0.9170 


0.261 


2 





250 


1.66 


23.56 


2.42 


(0.02, 7.93) 


1.1421 


1.1421 


0.265 


2 


+1 





0.00 


6.00 


7.35 


(6.46, 12.38) 


0.6569 


0.6569 


0.225 


2 


+1 


50 


0.58 


6.02 


7.28 


(6.41, 12.36) 


0.8484 


0.6587 


0.225 


2 


+1 


150 


1.52 


8.29 


5.66 


(5.11, 12.09) 


1.3162 


0.7473 


0.244 


2 


+1 


250 


2.21 


11.69 


4.29 


(3.50, 11.33) 


1.8062 


0.8580 


0.257 


2 


-1 





0.00 


6.00 


7.35 


(6.46, 12.38) 


0.6569 


0.6569 


0.225 


2 


-1 


50 


0.52 


7.46 


5.89 


(4.74, 11.58) 


0.5426 


0.7322 


0.239 


2 


-1 


150 


1.07 


19.16 


2.80 


( 0.68, 8.75) 


0.4937 


1.0626 


0.265 


2 


-1 


250 


1.37 


38.86 


1.66 


(-1.79,4.43) 


0.4315 


1.3797 


0.254 


2 


+2 





0.00 


6.00 


7.35 


(6.46, 12.38) 


0.6569 


0.6569 


0.225 


2 


+2 


50 


0.60 


5.18 


7.91 


(6.23, 11.88) 


1.0110 


0.6317 


0.213 


2 


+2 


150 


1.93 


4.52 


9.05 


(7.24, 12.22) 


1.7286 


0.5908 


0.202 


2 


+2 


250 


3.29 


4.31 


9.49 


(7.60, 12.32) 


2.4733 


0.5770 


0.198 


2 


-2 





0.00 


6.00 


7.35 


(6.46, 12.38) 


0.6569 


0.6569 


0.225 


2 


-2 


50 


0.51 


7.57 


5.81 


(4.63, 11.51) 


0.3583 


0.7375 


0.239 


2 


-2 


150 


1.21 


14.01 


3.59 


(2.08, 10.13) 


0.1992 


0.9384 


0.262 


2 


_2 


250 


1.66 


23.52 


2.43 


(0.04, 7.95) 


0.7552 


1.1411 


0.265 


r 


-1 


150 


2.61 


2.45 


16.67 


(11.81, 12.77) 


0.1336 


0.4352 


0.131 


r 


-1 


250 


3.20 


4.87 


9.00 


( 8.12, 12.92) 


0.3557 


0.5924 


0.208 



r -2 150 3.39 1.36 28.11 (14.05,11.73) 0.8025 0.3352 0.040 
r -2 250 3.90 3.27 13.37 (11.31,13.39) 1.4103 0.4859 0.164 
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The depth-dependence of eigenfunctions is determined primarily by the product A.B5. In 
particular, the radial orders in this model are given by n w 4.23 VAxS. For the selected 
modes they are between 26 and 34. Nonadiabatic parameters, / and T|, depend on n but 
also on the pulsation frequency in the star reference system, V sta r- 

The Hough ©-functions for selected modes are shown in Figs. 3 and 4. We may see 
the well-known effect of equatorial amplitude confinement, which increases with V rot . 
The effect is present in all modes including those with m—l though X(s) is a decreasing 
function in these cases. As for the remaining Hough functions, which are important, 
the confinement is also present. © and © have the same symmetry about the equator 
and the symmetry of © is opposite. 



6 Visibility of Slow Modes 

For the modes listed in Table 1 , we calculated amplitudes of light variation in the U and 
V Geneva passbands, Ay and Ay, with Eqs. (15), as well as that of the radial velocity, 
Ay ra d, withEqs. (19), adopting an arbitrary normalization, £=0.01. Coefficients ay and 
a s occurring in Eq. (15), were interpolated from the line-blanketed models of stellar 
atmospheres (Castelli and Kurucz 2004). 
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Figs. 5-7 show the Ay and Ay a d in function of the aspect angle. Because of the 
arbitrary normalization, they cannot be regarded as reliable predictors of expected am- 
plitudes. We may expect that chances of excitation are related to r| but it does not 
determine e. Still plots like these provide important hints for interpretation of the rich 
oscillation spectra such as that of HD 163868 (Walker et al. 2005b). For the t = 1 
g-modes we see a simple patterns in the dependence of Ay on the rotation rate. The 
aspect-dependence is qualitatively the same as in the case of no rotation. Rotation gives 
rise to a departure from the dipolar angular dependence of ^ r and 8 fboi but the contri- 
bution of the surface distortion to light variation remains small even at V mt = 250 km/s. 
Thus, the A y(i) -dependence is a simple reflection of the @(9)-dependence shown in 
Fig. 3. The equatorial confinement leads to a reduction of amplitude upon averaging 
contributions from the whole hemisphere. 
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Fig. 5. Amplitudes of light in the V-band of Geneva photometry, Ay, and of the radial velocity, Ay ra( i, at 
indicated equatorial rotational velocities plotted as functions of the aspect angle, for the I = 1 g-modes. The 
amplitudes are calculated assuming £ = 0.01 (see Eq. 1). 



The pattern of the dependence of radial velocity amplitude on the rotation rate, 
shown in the right panels of Fig. 5, is more complicated. In all three cases, the dominant 
contribution arises from pulsational velocity (the C pu i s term in Eq. 19). A secondary but 
significant contribution, the C m t term, adds in the case of prograde modes and subtracts 
in the case of retrograde modes. However, the main reason for the large difference 
between the modes and for the nonmonotonic aspect dependence is connected with 
properties of C pu i s - At high rotation rates the contribution from the advective term in 8v 
is large and causes that all its three components play a role. The large values of Ay ra( j in 
the case of £=l,m=+l mode for the near equatorial observers arise primarily from this 
term. We see that at least for some aspects, fast rotation increases the chances for mode 
detection by means of radial velocity measurements. The A y{i) -dependencies for the 
1 = 2 g-modes depicted in the left panels of Fig. 6, show that, despite of the increasing 
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equatorial confinement, the aspect of the best visibility moves toward the pole with 
increasing rotation rate. This somewhat unexpected feature, which has been already 
noted by Townsend (2003b), is explained in part by the effect of the surface distortion. 
There are aspects at which we see amplitude increase with V mt but on average the effect 
is the same as for £= 1, that is fast rotation decreases chances for photometric detection 
of slow modes. Also on average, amplitudes of the I — 2 modes are lower than those 
of I = 1. In the right panels of Fig. 6, we may see that radial velocity amplitudes of 
prograde modes increase with the rotation rate and that the effect is opposite for the 
retrograde modes. The role of the dot term is much more significant than at £ = 1. 
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We may see in Fig. 3 that ©(8) for the m = — 1 r-mode and for the £ = 1, m = g- 
mode look very similar. Likewise, in Fig. 4 we see the similarity of ©'s for the m = —2 
r-mode and for the £~2,m~ + l g-mode. Yet, as we may see in Fig. 7, the photometric 
amplitudes of the r-modes are much smaller and have different aspect dependence. The 
amplitude reduction is in part caused by cancellations in the integral over (j). Moreover, 
the effect of distortion, which is significant, cancels a part of the effect of temperature 
perturbation. The r-modes are antisymmetric with respect to the equator and are best 
seen from the intermediate aspect angles. The radial velocity amplitudes, shown in the 
right panels of Fig. 7, are much less reduced. Thus, spectroscopy gives a better chance 
for detecting r-modes. 



7 Prospects for Mode Discrimination 

Rotation has a very profound effect on slow mode visibility and hence on the procedure 
of mode identification. Unlike in non-rotating stars, the amplitude ratios depend on the 
azimuthal order, m, and on the aspect, i. Modes with various m differ not only in the 
surface geometry but also in their nonadiabatic properties. There are constraints on m 
following from localization in the frequency spectrum. Modes differing in the sign of 
m are well separated. As we have seen in Figs. 5-7, the aspect is an important factor 
in mode selection which has to be taken into account considering possible m and ; 
values. Nonetheless, a unique discrimination of modes is not likely possible without 
employing amplitude and phase data. 

Let us begin with photometric data alone. The observables, which do not depend on 
e, are amplitude ratios and phase differences from multiband photometry. In this case, 
the potential for mode discrimination rests mainly on the difference in the relative tem- 
perature and distortion contribution to the total light variation in different passbands. 
The latter contribution for the I = 1 g-modes is small, hence a discrimination between 
such modes might be difficult. However, as the plots in Fig. 8 show, these modes may 
be rather easily distinguished from the i = 2 and r-modes even if information about i 
and V rot is imprecise. A discrimination between the two I = 2 modes could be difficult 
if the measured phase difference is between 0.2 and 0.3. The range is consistent with 
m= — 1, if Vrot is between 150 km/s and 175 km/s, and m = — 2, if V rot is 250 km/s. 
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Fig. 8. Photometric diagnostic diagrams for selected modes based on light amplitudes and phases in the 
U and V Geneva passbands, for retrograde g-modes and for the r-mode with m = — 1 . The arrows indicate 
direction of increasing i. The i-ranges are limited by the condition A v > 10 mmag 



Any discrimination between these four modes based on frequencies alone would 
not be possible because all of them are retrograde and occupy the low frequency end 
of the oscillation spectrum. It is important to remember that the plots refer only to 
most unstable modes of each type. Thus, in real life, there is additional uncertainties in 
frequencies and in the /-values that affect the amplitude ratios and phase differences. 
With mean radial velocity data, discrimination between the two 1 = 2 modes should 
be unambiguous, as the plots in the lower panels of Fig. 9 show. In the upper panels 
we see that also distinguishing the two £ = 1 modes should be easy. Disentangling the 
m = +l and m = cases, however, may require data on line profile changes. 

The main benefit from radial velocity data is the possibility of simultaneous de- 
termination of I with no need for specifying the value of the complex parameter /. 
Instead, as described in the case of non-rotating stars by Daszyhska-Daszkiewicz et 
al. (2005), the combined spectroscopy and photometry data on amplitudes and phases 
may be used to determine /, which becomes an independent seismic observable, and 
the mode degree, £. In order to see how the method may be extended to the present 
case, let us rewrite Eqs. (15) and (19) as expressions for the complex amplitudes of 
light in the x-band and of radial velocity, respectively. From Eq. (15) we obtain 

* x (i) = <D'(i)fe + £'(i)e (20) 

where 

©'(/) = -1.086 ^Bi + aa), 

£*(;) = -1.086[2«i + « 2 - (2 + ar 1 )(a£Si + « 4 )] 

and e = eexp(imcpo). Similarly, from Eq. (19) the first moment of spectral line is equal 
to 

Mf (i) = i(S)R[9{ x {i)fe+ g x (i)e] (21) 
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Fig. 9. Diagnostic diagrams for selected modes based on amplitudes and phases in the V Geneva photometric 
passband and in the radial velocity, for g-modes with l= \ and 2. 



where 

r , . i sin/ ( \ r \ 

and 

§- x (0 = « 5 -|«7 +C5 ^ 6 -|« 8 ) -^5.1 [225 9 + «i0-(2 + C3- 1 ) (a^ 9 + «i 2 )] . 

These two equations are counterparts of Eqs. (1) and (2) of Daszyhska-Daszkiewicz et 
al. (2005). 

In the case of negligible rotation, the dependence of amplitudes on m and i has 
been absorbed in e and the problem was reduced to finding ef and e by the least square 
minimization of % 2 (l). The method requires stellar atmosphere parameters, which in 
fact may be improved on the process of pulsation amplitude fitting. An unique deter- 
mination of £ is possible even with imprecise atmosphere parameters if the min(x 2 )(^) 
dependence is strong. In the present case we have two more stellar parameters to im- 
prove which are i and V vo t- Moreover, we have to take into account the m-dependence. 
Prospects for mode identification depend on the strength of the the min(% 2 )(£,m) de- 
pendence. The plots in Figs. 8 and 9 suggest that it is likely the case. However, it 
remains to be seen when the method is applied to real data. 



8 Summary and Conclusions 

Our goal was to examine chances for detection and identification of slow oscillation 
modes whose frequencies are of the order of angular velocity of rotation. In our cal- 
culations of expected mode amplitudes, we relied on a nonadiabatic generalization of 
the traditional approximation, similar to that introduced by Townsend (2005a). The 
chances for detecting a particular mode depend, in part, on its intrinsic amplitudes, 
which may be calculated only in the framework of a nonlinear modeling. Our linear 
models give us only a hint which is the growth rate. Such models are expected to 
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be adequate for describing geometry of the mode which has important impact on mode 
visibility, its aspect and the angular degree dependence. With the formalism outlined in 
Sections 3 and 4, we calculated the observable amplitudes for selected unstable modes 
in models of a 6 M main sequence star as a function of the rotation rate and the aspect 
angle. The model may be regarded as representative for SPB variables. 

Departure from the individual spherical-harmonic dependence which increases with 
the rotation rate leads, in most cases, to lower photometric amplitudes. In contrast, the 
effect on radial velocity amplitude is most often opposite. However, the light to radial 
velocity amplitude ratio changes significantly from mode to mode and depends on the 
aspect. We showed, in particular, that the mixed r/g-modes are most easily detectable in 
radial velocity. Considering possible identification of peaks in rich oscillation spectra 
such as of HD 163868 (Walker et al. 2005b) it is important to take into account the as- 
pect dependence of considered mode amplitudes as it has been done by Dziembowski 
et al. (2007). 

The observables yielding numerical constraints on mode geometry and the aspect 
angle are amplitude ratios and phase differences. We calculated examples of diagrams 
based on photometric data in two passbands and on mean radial velocity measurements. 
Not always photometric data are sufficient for mode discrimination. Radial velocity 
data not only help discrimination but also allow to proceed in a less model-dependent 
manner. The photospheric value of the complex radial eigenfunction corresponding 
to the radiative flux may then be determined from data rather than taken from linear 
nonadiabatic calculations. A comparison of calculated and deduced values yields a 
constraint on the model. 
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Appendix 



Al. The Limb-Darkening Law 



For calculating the surface integrals in Eq. (15) and Eq. (19), we need to specify 
the limb-darkening law. Here, we use the nonlinear Claret (2000) formulae which we 
rewrite in the following form 



hx(fi) = 2 



1-2:4(1-^/2) 

k=l 

4 k 

1- 1 



The derivative with respect to jj can be easily obtained from this formula. The deriva- 
tives with respect to the effective temperature and gravity are given by: 
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A2. Integrals in the Expression for the Light Variation 

«i= J ©i'lsinGde 

o 

1+i 



where 



«2= J ^^sine:r : + /)i©'/> 3 j d9 
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o 

The function © is one of three Hough functions (see Section 2). The components of 
the unit vector directed to the observer (o r ,OQ,o^) are given in Eq. (13). 

A3. Integrals in the Expression for the Radial Velocity Variation 

2+' 

#5= J ©nsinGde 
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$6= J (©2'7 + ©2'8)de 
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